Cavity polaritons in the presence of symmetry-breaking disorder: closed-path time 

formalism 
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According to the mean-field theory of Zittartz, when subject to a symmetry-breaking disorder, 
the order parameter and the energy gap of an excitonic insulator are gradually suppressed up 
to a critical disorder strength. Recently, Marchetti, Simons, and Littlewood have used a replica 
trick to investigate the effects of disorder on the condensation of cavity polaritons. Within their 
nonlinear sigma model, it was found that the saddle-point equations assume the form reported 
previously by Zittartz in the contest of the symmetry broken excitonic insulator, but with an order 
parameter, to which both photons and excitons contribute. In this paper, we apply a closed- 
path time Green's function approach as an alternative to the replica technique to formulate a 
nonperturbative description of cavity polaritons in the presence of a symmetry-breaking disorder. 
A field theoretical method is used to derive the Schwinger-Dyson equations for the average photon 
field and the average single-particle Green's function. In contrast with the nonlinear sigma model 
and the corresponding saddle-point equations, we obtain that the exact Schwinger-Dyson equations 
cannot be mapped to the corresponding equations derived by Zittartz. This result not only shows 
that the theory of Zittartz cannot be applied to the excitons in a disordered quantum well coupled 
to the cavity photons with only minor modifications, but arises a question about the validity of the 
replica trick as well. 

PACS numbers: 71.35.-y 



I. INTRODUCTION 

The phenomenon of the Bose-Einstein condensation 
(BEC) in atomic gases and superconductors attracts 
much attention in recent years. Substantial efforts has 
been recently devoted to BEC of excitons in quantum 
wells (QW) and microcavities (MC)i. QW embedded 
within semiconductor MC have attracted considerable 
interest^ due to the following two reasons. Firstly, the 
recent progress in the growth and manipulation tech- 
niques of semiconductor heterostructures allows us to 
control the coupling between photons and excitons. The 
QW excitons embedded in semiconductor MC may be 
found in either weak- or strong-coupling regimes. In what 
follows we assume the strong coupling regime, where the 
photon-exciton interaction is larger than the exciton and 
photon damping rates, and therefore, the normal modes 
are mixed exciton-photon modes, called cavity polari- 
tons. Secondly, it is expected that the cavity polaritons 
should have bosonic behavior, and so are candidates for 
Bose condensation. 

Despite the progress made in semiconductor technology, 
a weak disorder may exist due to the following reasons: 
the interface roughness, the local thickness fluctuations 
during crystal-growth processes, randomly distributed 
impurities, boundary irregularities, and fluctuations of 
the alloy concentration of the epitaxial layers. The dis- 
order causes additional problems in the theory because 
all quantities of interest depend on the corresponding 
random potential, created by the disorder. The actual 
potential is unknown, but it is not important for the 
physical properties. Instead, the disorder is considered 
by means of the probability distribution of the random 
potential, i.e. one should perform the averaging over all 



possible random potentials. 

Turning our attention to the theoretical situation, we find 
that some authors^ have focused on the model, which as- 
sumes that the excitons are localized by the disorder, 
and can be described as two-level oscillators coupled to 
the light. In this model the original electron-hole-photon 
Hamiltonian is reduced to the Hamiltonian, which de- 
scribes the so-called generalized Dicke models The model 
is valid only for describing the very low energy excitonic 
states, and cannot be applied to a symmetry-breaking 
disorder, because it assumes that there exists only a 
symmetry preserving disorder potential, which is strong 
enough to localize the excitons. 

More complicated approach to the problem of excitons 
(or excitonic polaritons) in weakly disordered semicon- 
ductors is based on the assumption that the disorder af- 
fects only the center-of-mass motion, but does not af- 
fect the exciton internal degrees of freedom (see, e.g., 
Ref. [5[ and references therein). Recently, this idea 
has been applied to the cavity polaritons^. According to 
this disorder-affected-center-of-mass-motion (DACMM) 
approach, the two-particle Schrodinger equation for an 
isolated exciton separates into two equations: the Wan- 
nier equation for the relative motion of the electron- 
hole pair and the Schrodinger equation for the exciton 
center-of-mass motion in a random potential. As a re- 
sult, the two-particle exciton wave function can be factor- 
ized, and therefore, the coupling strength of an exciton to 
light is a random quantity, which depends on the exciton 
center-of-mass eigenfunctions ^(R). The next step in 
the DACMM approach is to use numerical simulations 
to generate random potentials. Once the energies £j 
and eigenfunctions ^(R) are calculated for a particular 
random potential, the radiative decay rates, the absorp- 
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tion (the optical density), or the exciton-photon coupling 
strengths can easily be evaluated numerically on a grid of 
a given number of points. There exists a many-body ver- 
sion of DACMM approach)! where instead of numerical 
simulations, the Green's function of the center of mass of 
an isolated exciton in the random field is calculated in the 
coherent potential approximation. It is expected that the 
factorization of the wave function is justified in the very 
low density regime. Strictly speaking, the DACMM ap- 
proach is based on the assumption that the disorder and 
the interactions (Coulomb and electron-photon interac- 
tions) could be treated independently. In other words, 
the factorization assumption greatly simplifies the prob- 
lem, but it separates the disorder and the interactions, 
and therefore, we may expect that the DACMM approach 
underestimates the influence of the disorder and may lead 
to incorrect conclusions. 

Decades ago, Zittartz^ demonstrated that the disorder 
and the interactions can be treated simultaneously if we 
perform the averaging over the disorder in the beginning 
of all calculations. Within this approach, the Green's 

functions are defined as < T {....} > and the brackets 
< ... > denote a thermal average, while / means the aver- 
age of / over the random potential created by the disor- 
der. Zittartz applied the Abrikosov and Gor'kov theory- 
developed for the case of superconductors in the presence 
of a symmetry-breaking disordered potential, to the case 
of an excitonic insulator in the presence of normal impuri- 
ties. The theory of Zittartz can be used with only minor 
modifications to investigate the effects of a symmetry- 
breaking disorder potential on the two-dimensional ex- 
citonic condensate in a high density regime, where the 
screened Coulomb interaction could be replaced by a con- 
tact interaction with a coupling strength g c . Assuming 
a Gaussian disorder potential with zero mean, and vari- 
ance V(r)V(r') = Ai5(r — r'), one can obtain the following 
set of equations for the order parameter A at the Fermi 
surface: 
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Here a = 2Am exc / A , m exc is the exciton reduced mass, 
u) m = (2m + 1)tt//3, (5 = (fcT)~\ where T and k are the 
temperature and the Boltzmann constant, respectively. 
The Matsubara summation in (Q]) must be cutoff at the 
energy which depends on the chemical potential fi. 
The solution of the above equations shows that: (i) the 
order parameter and the energy gap are gradually sup- 
pressed up to a critical disorder strength; (ii) the sup- 
pression of the energy gap is more rapid than that of 
the order parameter, which means that the existence of 
a gapless condensed phase is possible. 
Generally speaking, the case of cavity polaritons in the 
presence of a symmetry-breaking disorder is more compli- 
cated than the excitonic condensate, because the theory 



has to take into account the photonic contributions to 
all quantities of interest. To the best of our knowledge, 
there exists only one paper by Marchetti, Simons, and 
LittlewoocU^ (MSL), where a model for cavity polaritons 
in the high density regime in the presence of a symmetry- 
breaking disordered potential is proposed, treating the 
disorder and the interactions simultaneously. MSL have 
used the so-called replica trick to perform the averaging 
over the disorder. The replica trickii is based on the 

following relationship: InZ = limjv^o (Z N — 1)/N , 
where Z is the generating functional. Once replicated, 
MSL have decoupled the arising quartic term in Z N 
by means of the Hubbard-Stratonovich transformation 
with the introduction of a matrix field Q(t, uJ m ), and 
then, integrating over the fermionic fields the problem is 
reduced to the so-called nonlinear sigma-model action, 
previously used to study superconductors with magnetic 
impurities?^ The final step in this approach is to draw 
conclusions by investigating the structure of the saddle- 
point solution. At the level of the saddle-point approxi- 
mation the following three statements take place: 
(i) while the chemical potential does not exceed the cavity 
edge mode u c , an order parameter |A| ^ is developed: 
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(2) 



Here, g is the exciton-photon coupling constant, ip is the 
average photonic field, g\tp\ and |S| are the photonic and 
the excitonic contributions to the order parameter, re- 
spectively; 

(ii) the excitonic order parameter |E| and the photonic 
field | ip\ are not independent quantities because of the 
following constrain: 



(lu c - fi)\ip\ = (g/g c )\T,\; 



(3) 



(iii) the saddle-point equations in the high density regime 
can be mapped to the corresponding set of equations by 
Zittartz |T]), but with the order parameter |A|, defined 
by @, and with g c replaced by g eff = g c + g 2 /{u c - A 1 )- 
Because of the correspondence between the Zittartz 's 
equations and the saddle-point equations, MSL have 
concluded that in the low-density regime, where the ex- 
citations are mainly excitonic like (only a small fraction 
of photons contributes to the condensate), the order 
parameter and the energy gap are gradually suppressed 
up to a critical strength of the disorder. The suppression 
of the energy gap is more rapid than that of the order 
parameter, and therefore, the existence of a gapless con- 
densate is possible. When the density of the excitations 
is increased, the chemical potential rises (first linearly 
with the density), and when the chemical potential 
approaches u) c , the excitations become photonic like. In 
other words, the character of the condensate changes 
from being excitonic to photonic. 

The purpose of this paper is to show that all of the 
above conclusions are drawn only because MSL have 
performed the averaging over the disorder using the 
replica trick. To justify our point, we shall treat the 
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effects of a symmetry-breaking disorder on the cavity 
polaritons by applying the closed-path time (CPT) 
(or Keldysh) Green's function technique 1 ^. The main 
reason for using this approach is that in the case of a 
static random potential the Keldysh closed contour in 
the time direction leads to an automatically disorder 
independent generating functional. In other words, the 
CPT approach allows us not only to avoid the need to 
introduce replicas, but to perform the averaging over 
the disorder in the beginning of all calculations as well, 
which is the main requirement when the disorder and 
the interactions are treated simultaneously. The special 
form of the Keldysh time contour automatically ensures 
that the denominator in the representation of the Green 
functions via functional integrals is equal to unity. The 
last allows us to derive the exact equations for the 
average photon field and for the average single-electron 
Green's function. In the quantum-field theory, these 
equations are known as the Schwinger-Dyson (SD) 
equations. In contrast with the saddle-point equation 
(J2|) , the exact SD equations clearly indicate that in the 
presence of a symmetry-breaking disorder the photonic 
contribution to the mass operator is not proportional 
to the average photon field. We shall see that the 
results derived by applying the replica trick correspond 
to the assumption that one can replace the average of 
the product of two random functions with the product 
of the corresponding average functions. However, it 
is known that the average of the product is not the 
product of the averages, and therefore, the validity of the 
conclusions based on the replica trick is questionable. 
The exact result that the photonic contribution to the 
mass operator is not proportional to the average photon 
field does not allow us: (i) to map the SD equations 
to the corresponding Zittartz equations; (ii) to use 
the Ward identities in order to prove the existence of 
the Goldstone mode below the critical temperature, 
and therefore, the question about the existence of a 
condensate of cavity polaritons in the presence of a 
symmetry-breaking disorder remains open. 
The remainder of the paper is organized as follows. In 
Sec. II, we discuss the formation of a condensate in 
the absence of a disorder. This is because we intend 
to check the validity of the saddle-point approximation 
by eliminating the effects generated by the replica 
trick. By applying the Matsubara Green's function 
method we demonstrate that the polariton spectra can 
be obtained from the common poles of the photon and 
the two-particle electron-hole Green's functions. It turns 
out that the saddle-point equations in the absence of a 
disorder lead to the same conclusions as those drawn 
by applying the Matsubara Green's function method. 
The approach used in Sec. II is very general and, in 
principle, it is able to treat any density regimes. It 
also allows us to prove the existence of the Goldstone 
mode below the critical temperature. Our method 
provides a set of coupled BCS and Bethe-Salpeter (BS) 
equations similar to the corresponding equations for 



an excitonic condensation ^ 4 ' 15 Due to the photonic 
contribution to the condensate, the BCS and the BS 
equations are more complicated than the equations 
reported in our previous paper—, and it would be a 
very challenging task to solve them in the case of a 
low-density limit. Such an ambitious task will be left 
as a subject of future research. In Sec. Ill we treat the 
effects of a symmetry-breaking disorder on the cavity 
polaritons by applying CPT Green's function technique, 
because this approach is analytical nonpertubative one 
which provides exact results. Furthermore, the Keldysh 
formalism could be applied to the cavity polaritons in 
noncquilibrium conditions. 



II. CAVITY POLARITONS IN THE ABSENCE 
OF A DISORDER - THE MATSUBARA GREEN'S 
FUNCTION APPROACH 

The system under consideration consists of a single 
QW grown inside a semiconductor MC is an arrange- 
ment of two-plane parallel mirrors with reflectivity close 
to unity. The two infinite and parallel perfect mir- 
rors are perpendicular to z-axis, separated by a dis- 
tance Lq, one mirror is at z — Lq/2, and the other at 
z = —Lq/2. In what follows we are interested in the 
case of a single QW extending over — L/2 < z < L/2 
made from a direct-gap semiconductor with nondegener- 
ate and isotropic bands when the electron-hole motion 
along the z-direction is confined between two parallel, 
infinitely high potential barriers. With the perfect con- 
finement approximation the dispersion laws for electrons 
and holes are E c (k c , A) = E g + k 2 /2m c + tt 2 A 2 /2m c L 2 
and E v {k v ,^) = —k 2 v /2m v — 7r 2 £ 2 /2m v L 2 , respectively. 
Here m c (m v ) is the electron (hole) effective mass, E g is 
the energy gap, and k CilI is a two-dimensional (2D) wave 
vector. A, £ = 1, 2, ... denote the quantum number of the 
states in the infinitely deep wells. In what follows we use 
the simplest approximation which takes into account only 
the first electron and hole confined levels, i.e. A = £ = 1. 
For each photon wave vector there are two possible po- 
larizations: one with transverse electric field (TE), and 
second, with transverse magnetic field (TM) In what 
follows we will take into account only the TE modes 
which interact with transverse polarized excitons. The 
longitudinal photon modes mediate the Coulomb inter- 
action between the charges in the QW, but we neglect 
this effect assuming that the Coulomb interaction be- 
tween the charges in the QW is affected only by the con- 
finement of the charges. The cavity- mode dispersion is 
^s(q) = cy/q 2 + (tts/Lq) 2 , where s = 1,2,..., and q is 
a 2D vector. In the following, we suppose that the only 
s = 1 cavity modes f2(q) interact with the electron sys- 
tem. 

In terms of the field theory, the transverse and the lon- 
gitudinal photon modes are described by boson fields 
^4_l(p) and A||(p), respectively. They interact with the 
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electron system, described by fermion fields ip + (y) and 
ip(x). The total action of the system is 

s = s { e) + S^ ] + s^l 

The actions for non-interacting electrons and photons are 



and 



(0)-l 



+ ±A ± (p)D^- 1 (p lP >)A ± (p>), 



respectively, 
scribed by 



The electron-photon interaction is de- 



i/>(y)r\°\v,x\pMx)A\\(p) 

+ ^(y)T^(y,x \ p)i,(x)A x (p). 

The composite variables y = {r, u}, x = {r',u'}, and 
p = {R, v} are defined as follows: r, r',R are 2D ra- 
dius vectors, and according to imaginary-time (Mat- 
subara) formalism the variable u,u',v range from to 
h[3 = h/(kT). We set h = 1 and we use the summation- 
integration convention: that repeated variables are 
summed up or integrated over. G^- 1 (y,x) is the in- 
verse single-particle Green function for non-interacting 
electrons in a periodic lattice potential G^~ 1 (y, x) = 
E Um e ""'" ,( °"°' ,G(0) " 1 ( r , z 1 r ' z '; ! w m )- The function 
G , ^°^~ 1 (r, z, r', z 1 ; iuj m ) is defined as a sum of elec- 
tron J2u c y*,k c ( r ' z )¥ , c,k c (i" / ,z') G '£'' _:L (k c ;icj m ) and hole 
Ek„ ^*.k„( r ' z )^«,k„(r',^')G'l° ) ~ 1 (k 1 ,;iw m ) parts. Here 
Gic )_1 (k c ; iu) m ) = iu m - [£? c (k c , A = 1) - Li c ], and 
G^J^ 1 (k v ;KJ m ) = iuj m - [E v (k v , £ = 1) - Hi,}. The func- 
tions iy5 Cj k c ( r ) and <p Vi \s, v (r) are the wave functions of 
the first electron and hole confined levels, defined by 
the solutions of the corresponding Schrodinger equations. 
The electron and hole chemical potentials are denoted 
by \i c and p v , respectively, and the symbol XL m i s 
used to denote f3~ 1 J2 m - For f erm i° n fields we have 
uj m = {2n/f3)(m + 1/2); m = 0, ±1, ±2, .... 
In addition to the lattice potential, the electrons and 
holes experience a Coulomb interaction, described by the 
termrjj >Z?j >rj°J: 



r^(y 2 ,x 1 \p)D\;'>(p,p')r^(y 3 ,x i \p') 



,(0) 



(0), 



S(m - u 3 )6(u 2 - Ui) V y~] ipi^ z {Yi)Lp* M ir 2 )x 



E ^ 



Here i,j — {c, v}, is the longitudinal part of the 

photon propagator (in a gauge, when the scalar po- 
tential equals zero) and is the vertex. Vb(q) = 



27re 2 /(X|q|)/|q| denotes the Fourier transform of the 2D 
unscreened Coulomb potential. The structure factor /(x) 
takes into account the first confined QW electron and 
hole levels: 



/(*) 



3x 2 +8ir 2 32tt 4 [1 - exp(-x) 



x{x 2 + 47T 2 ) x 2 (x 2 + A-k 2 ) 2 

The inverse transverse photon propagator is: 



,t[q.(R-R')-u p (v-v')] n(°) _1 



D^- 1 {^uj p ). (4) 



q u p 



Here A Q is the area of the cavity, and 1 (q,iu> p ) — 

2ttc 2 /{(iLUp) 2 — f2 2 (q)]. The symbol J2 Up is used to denote 

/3~ 1 J2 P - For boson fields we have uj p — (2ir//3)p;p = 
0, ±1,±2, ... We assume that the electron-hole-photon 
system is in thermal equilibrium, which means that the 
poles of the photon retarded Green's function have to 
be obtained from the corresponding Matsubara Green's 
function, by the substitution iuj p — ► u> + p, + i0 + . Here 
fx = ji c — ji v is the chemical potential of the system. The 
vertex (y, x | p) has the following form: 

T { °\y 2 ,xi | p) = Tf ) (r 2 ,u 2 ,ri,ui K,v) = 



8(ui - v)5(ui - u 2 ) 



E E ^vi^vi.M 

q i.kij.kj 



< Jjkj I j(q)-n(q) | i,kj >, 



(5) 



where j(q) denotes the single-particle current operator, 
and n(q) = e q x e z , where e q = q/q and e z = (0,0, 1). 
It is well-known that all Green's functions can be ob- 
tained by functional differentiation from the generat- 
ing functional W[J, M] = \nZ[J,M] of the connected 
Green's function, where 



Z[J,M] = J Dii®,il>,A ± ,A\ { ]x 

exp[S + J ± {p)A ± (p) + J,, (p)A l{ (p) - ^{y)M(y, x)^{x)] . 

(6) 

Here J and M are the sources of the corresponding fields. 
By means of the functional (|6|) we introduce the following 
functions (after the functional differentiation one should 
set J=M=0): 
Average photon field: 



R\\,(A-)(P) 



sw 



^II(-L)O) 



single-particle Green function: 



G(x,y) 



6W[J, M] 
'SM(y,x)' 



(7) 



(8) 
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transverse (longitudinal) photon Green function: 
5 2 W[J, M] 



*%x(p). 

A ' W 



5J ±( ||)(P')^±(||)(P) <%x(/>') 
two-particle electron-hole Green function: 



K 



x y 



8 2 W[J,M] 



SG(x,y) 

' 8M(y',x')5M(y,x) 8M(y',x')' 



(10) 

transverse (longitudinal) electron-photon vertex func- 
tion: 

rx (B )(y,*|p) = -^^^„,a/,P). (ID 

As a consequence of the fact that the measure is invariant 
under the translations ip — > i/j + Si/j, Ax,« — > -^X,ii + <^4x,u 
we derive the SD equations 



J\\(±)(p) ~ ^(X) 1 ^' P') R M±)(P') 
+ rf^ ) (y,x\p)G(x,y)=0, 



(12) 



G~ l (y,x) - G^- 1 (y,x) + M(y,x) + E(y,x) = 0. (13) 
The mass operator £ has the form: 

£(y, x) = (y, x\p)R± (p) + rj 0) (y, x \ p)R {l (p) 
-r^(y,x'\p)G(x',y')r ll (y',x\p')D ll (p,p') 
-r<£\ y ,x' | p)G(x',y')T x (y',x \ p')D x (p,p'). (14) 

The charge neutrality leads to the photon held J?m (p) 
that vanishes identically. Thus, all terms proportional to 

R\\{p) (or D\\(p,p')T^\y,x \ p')G{x,y)) does not need 
to be taken into account because of the global neutrality 
of the electron-hole system. In what follows we assume 
the so-called Hartree-Fock approximation in which the 
mass-operator has the form 



^(y,x) = T { l\y,x | p)R±(p) 



r\°\y,x> | p)G(x',y')T ll (y',x \ p')D^p,p'). (15) 



The first and the second terms in (|T5|) are called 
the Hartree term and the Fock term, respec- 
tively. The presence of Bose-condensed polaritons 
modifies the single-particle Greens functions, and 
therefore one has to consider the so-called nor- 
mal G~J-(k;icj m ) = G£ )_1 (k;iw m ) - E cc (k) and 

G~^(k;tu> m ) — G™ 1 (k;tu> m ) — E^k), and anomalous 
G~^{k;iuj m ) = G~^(k;iuj m ) = -A(k) single-particle 
Green's functions (the spin degrees of freedom are not 
included). The diagonal parts of the mass operator in 
the Hartree-Fock approximation are as follows: E cc (k) = 
E q ^(k - q)E, m G cc (q;^ m ), E w (k) = E q U(k- 
<l)ZL G vv (q;iuj m ). Here U(q) = 27re 2 /(L|q|)/e 00 |q| 



denotes the screened Coulomb potential. Using (| 1 5|) we 
calculate for the non-diagonal parts of the mass operator 
in the Hartree-Fock approximation: 

A(k) = ^r ( L 0) (q,k) J R x (q) + A exc (k), (16) 
q 

where A is the order parameter for the system. The first 
(Hartree) and the second (Fock) terms in l|16p represent 
the photonic and the excitonic contributions to the order 
parameter, respectively. The exact form of V ® (q, k) can 
be calculated by means of ^ and ©, but we assume 
that the photons are coupled to the electron-hole system 
through the local interaction, i.e. r^(q, k) = g(q — 
k). In this approximation the photon field R± and the 
excitonic order A exc parameter are defined as follows: 

R ± (k) = S(k - q) Yl G -(q; w m )/ (0 2 (k) - p 2 ) , 

(17) 

A exc {k) = v(k - q) £ G -(q; «"m)- (is) 

The photonic and the excitonic order parameters are not 
independent. Using Eq. (|18[) for the excitonic order pa- 
rameter, we calculate J2u m G C v((i; «w m ) = '%2 k V~ 1 (ci — 
k)A exc (k), and therefore, we obtain the following con- 
straint: 

*J-(k) = n>ul ,. 2 Eg( k ~ ^"'fa - P)^e XC (p), 



n(ky 



q.p 



where we have introduced a function V 1 defined by 
E p U(k-p)U- 1 (p-q) = <5(k-q). 
In the Hartree-Fock approximation the normal single- 
particle Green's functions are>^ 



G cc (k; iuj m ) 



G vv (k;iLo m )) 



w+(k) 

„,2 



■ «_(k) 



_ KV m - LU+ (k) 

G cl) (k; iuj m ) = G vc {k; iu) m ) 

I 1 



o;_(k) 



(20) 



_ iu m - lu + (k) iu> m - to- (k) 
Here, the following notations have been used: 



u 2 (k) 



1 



ry(k) 



s(k) 



,v 2 (k) 



1 



V (k) 



e(k) = V^ 2 (k) + A 2 (k) 



e(k)_ 

w±(k) = C(k)±e(k), 



C(k) = - [E c (k, 1) + ^(k, 1) - /i c - p v ] 
+ l^U(k-q) [n_(q)-n + (q)], 



77(k) 



[.E c (k,l)-.E B (k,l)-d 



-i]TV(k-q) 



l-[l-n+(q)-n_(q)] 



??(q) 



e(q) 
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where n± (k) = [1 + exp (±/3w±(k))] L . Below the 
critical temperature the order parameter is developed 
(A(k) ^ 0) and the single-particle excitations are co- 
herent combinations of electron- like w+(k) and hole- like 
cj_(k) excitations, renormalizcd due to the interaction 
with the cavity modes. The coefficients u(k) and i>(k) 
which are called coherent factors, give the probability 
amplitudes of these states in the actual mixture. 
It is expected that the BEC phenomenon is not sensitive 
to the difference in electron and hole effective masses, 
and therefore, we assume m c = m v = 2m exci where 
m exc = m c + m v 1 i s the exciton reduced mass. The 
equal mass assumption simplifies very much the calcu- 
lations because in this case p c + Pv = E g , ((k) = 0, 
n+(k) = n_(k) = [1 + exp (/3e(k))]~ , and rj(k) is de- 
fined by: 



T,(k) 



E„ 



2 ^ 
q 



V(k - q 



2mL 2 2m exc 



1 — tanh 



?7(q) 



£(q) 



(21) 



Note, that in our approach ry(k) is defined self- 
consistently by the solution of Eq. (|2~lj) . while MSL 
neglected the contributions due to the diagonal parts 
(S cc (k) and T, vv (k)) of the mass operator, using the fol- 
lowing expression: 



?7(k) = k 2 /4m exc - e F , 



(22) 



where the effective Fermi energy is defined by Ef — 
p 2 F /4m exc = (p - Eg - TT 2 /2m exc L 2 )/2. In the high- 
density limit, the neglected diagonal parts of the mass 
operator are responsible only for small renormalization 
of the single-particle excitations co±(k). But, in the low- 
density limit the diagonal parts of the mass operator are 
crucial, and therefore, they cannot be neglected. 
The order parameter ([To]) includes both excitonic and 
photonic contributions: 

A(k) = ]T 5 (k-q)i? ± (q) + A e;cc (k), 



and is determined by the constraint and the follow- 
ing BCS self-consistent equation for the excitonic order 
parameter: 



A e£CC (k) = ^(k-q)£ 



A(q) 

+ e 2 (q)' 



(23) 



(i2j_(k) = Rdk t o, and A exc (k) = A exc S kfi ), we obtain 
the total order parameter A = gR + A exc , where 



A exc — g c ^ ^ 



A 



k w„ 



e 2 (k)' 



While the chemical potential p does not exceed the 
cavity-mode energy uj c at k = 0, the constrain (JTUJ) as- 
sumes the form: 



R 



9 



9c (u> 2 -p 2 ) 



(24) 



Similar relationship has been found by MSL, but with 
uj c — p in place of ui 2 — p 2 . This is because of the dif- 
ferent photon Green's functions, used by MSL. In the 

2 



high-density regime 77 (k) = k /m exc — £p, and Eq. (j 231 
assumes the form of the BCS gap equation for a super- 
conductor: 



1 



9 



- P 2 



E 



1 



2e(k) 



tanh 



/fe(k) 



(25) 



When 3 = 0, the last equation assumes the form JT]) with 
a = 0. 

The chemical potential p is a nontrivial function of the 
total number of excitations N = N p h + N exc in the con- 
densate and should be calculated by solving the BCS 
equation self-consistently. The total number of photons 
is Nph = J2k ^xO 4 )- The number of condensed electron- 
hole pairs N exc is: 



N exc = J2 



1 — tanh 



(3e(k)\ r?(k) 



e(k) 



(26) 



Our next step is to show that the photon Green's function 
and the two-particle electron-hole Green's function have 
common poles - the excitation spectrum in the presence 
of a condensed phase. To prove the last we introduce a 
Legendre transform 

V[R,G] = W[J,M]+J a (p)R a (p)+M(y,x)G(x,y). (27) 

The repeated Greek index a denotes summation over the 
parallel || and the perpendicular _L components of the 
corresponding quantity. By means of the above Legendre 
transform, we derive the following exact equations^ 

D ± (p,f/) = Df{{p,p')+Df{p,p")Tf([y,x \ p") 



K\ * I )T^(y' 1 x'\p"')D^(p"',p'), 



(28) 



To make contact with the equations of MSL in the ab- 
sence of a disorder, we assume a high-density limit. In 
the high-density regime the screened Coulomb interac- 
tion could be replaced by a short-range contact interac- 
tion with a coupling strength g c given by angular average 
over the Fermi surface. Assuming that: (i) ,g(q) = .9, 
and (ii) the order parameters are space independent 



SG(x,y) 
SJ ± (p) 

Ki0) ( X y ^)r±(jAa/|p')£±(p',P)= (29) 
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y x 

x y' 



= K (o)-i y x )_j 

1 x y ' 



y x 
x y' 



(30) 



where f J J, j = G(x, y')G(x', y) is the free two- 
particle propagator and the kernel of the BS equation 
(|5D|) is given by: 



f y x'\ = 8Y,{y,x) 
\x y' ) 5M(y',x>) 



') (31) 
+ T^(y,x\p)D^(p 1 p')T^(y' 1 x'\p'). 

There are two important conclusions that could be ex- 
tracted from the above equations. The first one follows 
from Eq. (|29]) . This equation clearly indicates that the 
Fock term in Eq. (fT5|) is related to the two-particle 
Green's function, so we can write the mass operator in 
the following form: 



E(y,x) = T^(y,x\p)R a (p) 

-r£%,z'|p)tf(*; y J, 

rW(y" > x"\p')D(V(p,p')G- 1 (y',x). 



(32) 



The second conclusion follows from Eq. (|28)) . Obvi- 
ously, the transverse photon Green's function and the 
two-particle Green's function have common poles - cavity 
polaritons. The poles of the transverse photon Green's 
function are defined by the solutions of the Maxwell equa- 
tions for a transverse wave e(Q,u/) = £l 2 (Q)/o; 2 , where 
the dielectric function e(Q, oj) in the case of a single ex- 
citonic resonance at energy £;(Q) is: 



d 



Ef{q)-uj 2 - 2iw 7o 



Here t b is the background dielectric constant, 70 is the 
broadening of the excitonic resonance, and d is propor- 
tional to the corresponding oscillator strength. In prin- 
ciple, the excitonic resonance energy and the oscillator 
strength can be calculated by solving a set of coupled 
BS equations for the energy and the wavefunctions of 
the quantum- well excitons. These equations are similar 
to equations (42) and (43) from our earlier paper The 
only change that should be done is related to the exis- 
tence of an extra term, Ty DyTy G, in the mass oper- 
ator. The last term contributes to the exchange interac- 
tion in the similar manner as the term rj| rj| G gen- 
erates the analytical exchange interaction between elec- 
trons and holes. The exchange interactions are important 
only for the fine structure of exciton levels, and therefore, 
we neglect the exchange interaction terms. As a result, 
we obtain the BS equations similar to the case of the ex- 
citonic condensate. Generally speaking, to calculate the 
excitation spectrum in the presence of a condensed phase, 
one has to solve the BCS and the BS equations simulta- 
neously, taking into account the fact that the chemical 
potential depends on the number of excitons and photons 



in the condensed phase. Such an ambitious task will be 
left as a subject of future research. 

We finish this section with a brief discussion of the so- 
called Thouless criterion^ In the case of superconduc- 
tivity this criterion says that below the critical temper- 
ature T c the T-matrix has a pole at zero frequency and 
zero momentum (the existence of Goldstone mode be- 
low T c ). To check whether the Thouless criterion works 
in the case of condensed cavity polaritons we follow the 
method based on the Ward identities*^ Taking into ac- 
count the fact that in the case of MC polaritons the 
mass operator depends on the photon field and 
the Green's function G, we first invert the SD equations 
to express the sources J a (p) and M(y,x) as functionals 
of the field R a and the Green's function G. Second, we 
assume that there exists a continuous transformation, for 
example, a rotation in order-parameter space, which de- 
pends continuously on the parameter A. If the system is 
invariant under this transformation, then the variation 
of the Legendre transform implied by the transformation 
is equal to zero, i.e. 5\V = 0. The BS equation for 
the two-particle Green's function K = + K^°'IK 
can be rewritten in terms of the many-particle T-matrix, 
T = I + IK^T, in the form K = K^+K^TK^. By 
means of the last form of the BS equation and the defini- 
tion (fTI))) we calculate the variation of the inverse Green's 
function SxG- 1 = - (l + TK^) S X M . Using the SD 
equations we find 5\G~ l = —8\M — 5\E, and therefore, 
T _1 5aS = K^S X M. But, according to (J57J) we calcu- 
late 8\M = ^ (5 X V) = 0, and therefore, T^SxE = 0. 
Above the critical temperature T c the order parameter 
is zero, and hence, S\T, = 0. Thus, T~ 1 5\I] = is sat- 
isfied trivially. Below T c the order parameter is nonzero 
and 8\Y> ^ 0, which requires that the inverse T-matrix 
has a zero eigenvalue. Thus, we conclude that below the 
critical temperature T c , the T-matrix must have a pole 
at zero energy and zero momentum. The existence of 
Goldstone mode below T c indicates that the formation of 
a condensate in MC is possible. 



III. CAVITY POLARITONS IN THE 
PRESENCE OF A SYMMETRY-BREAKING 
DISORDER - THE KELDYSH GREEN'S 
FUNCTION APPROACH 

In the presence of a disorder we define the Green's 
functions as an average of the time-ordered products of 
the fields, but the average includes both the quantum 
and the disorder averaging. We use < F > for the quan- 
tum averaging, and F for the disorder averaging. The 
Matsubara Green's function approach discussed in the 
previous Section, cannot be directly applied to the cavity 
polaritons in the presence of a disorder, because one has 
to calculate < InF >. Nevertheless, for a given disorder 
configuration we can write the relationship between the 
random photon Green's function G and the random two- 
particle Green's function K, similar to Eq. 11280 . After 
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that, the disorder averaging replaces the random Green's 
functions on the both sides of Eq. by their aver- 

ages. Thus, we obtain that average photon D and aver- 
age two-particle Green's function K have common poles 
- the cavity polaritons in the presence of a disorder. The 
next question to be answered is about the possibility to 
observe a condensate in MC in the presence of a disor- 
der. To answer this question we have to examine the 
non-diagonal parts of the average single-particle Green's 
function. The last can be obtained by solving the SD 
equations. We shall use the Keldysh technique to derive 
the SD equations in the presence of a disorder, because 
the CPT approach allows us to perform the average over 
the random potential exactly, i.e. non-perturbatively. 
Within this approach we have boson (photon) longitu- 
dinal and transverse fields A a (z) = Am x(R, t^) inter- 
acting with a fcrmion (electron) field tp + {y) — i/jj~(r,t), 
or ip(x) = ipj'(r' ,t'), in the presence of disorder. The 
variables y,x and z are composite variables y — {r,i, j}, 
x = {r',t',j'}, z = {R, t"}, where r,r',R are the corre- 
sponding 2D radius vectors. The index j = 1 (or j = c) 
denotes the electron states, and j = 2 (or j = v) de- 
notes the hole states. For simplicity, we suppose that the 
electrons are spinless. At a zero temperature the total 
action of the system is S = + + To 
incorporate the effect of disorder we introduce random 
static symmetry-breaking (charge-dependent) potential 
V(r). In the presence of a disorder the action which cor- 
responds to the electron system assumes the form: 

= tp + (y)^ ^ 1 (y, x) - V(y,x)]iP(x), 

where V(y,x) = V(r, t, j; r',f ,/) = V(r)S(r -r')5(t- 
t')djji . The actual potential V(r) is unknown, but we 
assume that it obeys Gaussian statistics such that: 

V(fj=0, V(r)V(r<) = AS(r - r'), (33) 



where A = {2-kvt)^ 1 , v = m exc /ir is the 2D density of 
states, and r is the corresponding scattering time. The 
disorder averaging in (|33| is defined for any functional 
F[V] by the functional integral: 

F = J DVF[V]exp[~^ J V 2 (r)dr}. (34) 

The actions and S^ e ~^ are given by: 

= \A a (z)DW- 1 (z,z')A a (z>), 

S (e " w) = ^ + (y)r^(y,x | z)^{x)A a (z). 

In the CPT formalism, the Green's functions are defined 
by means of two time orderings in the same formula. In 
what follows we use the single-time representation. In 
this representation the two time orderings are replaced 
by a single time ordering along the Keldysh contour 
which at a zero temperature consists of two branches: 
the right-going (+) from — oo to oo and the left-going 
(— ) from oo to — oo. The symbol t means that the time 
integral J dt along the Keldysh contour could be written 
as two usual integrals, i.e. J dt = dt + — dt~ . In 
other words, the time variable t on the positive branch 
equals t = t + , and t = t~ on the negative branch. It 
should be mentioned that our equations are valid for 
both nonequilibrium and equilibrium conditions, but we 
do not discuss time-dependent phenomena on an ultra- 
fast scale. Instead, we intend to describe steady-state 
phenomenon, such as the light propagation in crystal in 
terms of excitonic polaritons. In the steady-state regime 
all quantities depend on the relative time t — t" . 
The inverse free propagator G^~ 1 (y, x) = 

Gjj~ 1 ( r ,i; r', t') is defined as follows: 



/OO 7 

K ~°° 

/OO J 



(35) 



k„ 

Here Gc°' )_1 (k c ; lu) = uj — [E c (k c , A = 1) — fi c ] + ^0 + , 1 (k v ;oj) = oj — [E v {k v , £ = 1) — /z„] — tO + are the inverse 

zero-temperature free electron and hole propagators. 

In addition to the lattice and the random potentials, the electrons and holes experience a Coulomb interaction, 
described by the term rf.^D^r^: 

r^^xill^^fezOr^^Xik') =<%-t3>%-*4)£ E ^o(q)^,k J (r 1 )^ kj+q (r 2 )x 

J-:' k P ' •! (36) 

<Pj', Vjl (i-3)^', Pj ,- q (r4). 

Here -Dy is the longitudinal part of the photon propagator (in a gauge, when the scalar potential equals zero) and 
r||°' ) is the vertex. Vb(q) denotes the Fourier transform of the 2D bare Coulomb potential, and has been defined in 
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Sec. II. 

The inverse transverse photon propagator is: 



D^-\z,z>) = Df - 1 (R ) i;R',f) = ^^E /°° f^^-^" 1 ^), (37) 

Here A is the area of the cavity, and £>^ _1 (q,cj) = 27rc 2 /[(u — ^) 2 — 2 (q) + «0 + ]. The vertex T±\y,x \ z) has the 
following form: 



rf (m,xi | *) = r^,(r 2 ,i2,r 1; ti | R,t) = 



_ r (0) 



E E eiq ' R ^,k 3 (i- 2 )^, Py (ri) < j.kj | j(q).n(q) | j', Pj , >, 



(38) 



q k ,p ■ 



r 



where j(q) denotes the single-particle current operator, 
and n(q) = e q x e z , where e q = q/q and e z = (0, 0, 1). 
Let us introduce the generating functional W[J, M; V] of 
the connected Green functions: 



W[J, M; V] = -ilnZ[J, M; V], 



(39) 



where J = J\\,±(z) and M = M(y,x) are the sources, 
and the functional Z[J, M ; has the form: 

Z[J,M;V] = /" D/iexp{i[S 

+ J\\{z)A\\(z) + J±(z)A ± (z) - ^ + (y)M(y,x)^(x)}}. 

(40) 

Here D/j, = CDtp + Dtp DA denotes the functional mea- 
sure. The success of the Keldysh technique is based on 
the fact that the normalization constant C is disorder- 
independent. Thus, we assume that C is chosen in the 
manner that Z[J = 0,M = 0;V] = 1. It is clear that 
because J and M do not have the same behavior on the 
forward and backward parts of the Keldysh contour, the 
generating functional is not equal to unity if the sources 
are not nullified. 

By means of the generating functional of the connected 
Green functions wc introduce the following average quan- 
tities: 

Photon field R a (z_) (in what follows a =||, or a =-L): 



SW[J,M;V] 



SJ a (z) 

photon Green's function D a (z,z'): 



S 2 W[J,M;V] 



5J±(\\)(z)SJ ±m (z>) 
single-particle Green's function G(x,y): 



J=M=0 



(41) 



(42) 



SW[J,M, V] 

G{x,y) = —i ira~r7~~~\ \j=m=o 



5M(y,x) 

/XV' 

two-particle electron-hole Green's function K [ ~ —, 

y 



(43) 



K 



x y 
y x' 



S-W .1. U:\" 

' 5M{y,x)6M(y',x'y J=M=0 ' 



(44) 



Evidently, we have four single-electron Green's functions 
Qmv2 anc j f our photon Green's functions D Vin2 , where 
r]i,r)2 = + or — depending on whether the time variable t 
is on the positive branch or on the negative branch. Note, 
that the time integration in expressions like C(x 1 ,x 3 ) = 
A(x l7 y 2 )B(y 2 ,x 3 ) follow the convention: 

C^{x 1 ,x 3 ) = C{v 1 ,tf;v s ,tf) 
= I dr 2 r dt+A(r u tf ;r 2 ,t+)B(r 2 ,t+;r 3 , tf) 

J J — OO 

dta A(n , tf ■ r 2 , q ) B (r 2 , q ;r 3 ,tf). 



dv 2 



Due to the 
K ( n,*i r 3 ,i 3 
r 2 ,i 2 r4,i 4 



(45) 



timc-translational 



depends on t 12 = — 1 2 , t i3 — t 4 



mvanance 



and i 31 = t 3 — t 1 . In what follows, we shall see that 
our equations will involve the two-particle Green's 
functions with t 12 



0, and therefore, we 



have four different two-particle Green's functions 



K 



r-utT r 3 , tj 



The corresponding re- 



r 2,ii r4,i 3 

tarded Green's functions G R ,D R ,K R , for example, G R , 
can be expressed as G R = G - G~+ = G+~ - G++. 
The Keldysh technique allows us to perform the disorder 
averaging. The resulting equations for the average pho- 
ton field R a (z) and the average single-particle Green's 
function G(x, y) are as follows: 



6Z[J,M] I 

RaU) - - l ^UzT lj=M = ' 



5Z[J,M] ^ 



(46) 



(47) 
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where the average generating functional Z[J, M] = Z[J, M; V] is defined by the equation: 
Z[J,M] = [ D f iexp{ l [i; + (y)G^- 1 (y,x)i:(x) + ^A^D^' 1 (z, £)A a {£) + ^ \y)Yf> (y,x\ z)^{x)A a (z) 



(48) 

+ ^AV^ + (r,i)^(r,t)^t(r,i')^.(r,i') + J a (z)A a (z) - 4> + (y)M(y,x)^(x)}}. 



To calculate R and G one has to know the functional Z[J, M}. Note that Z[J = 0, M = 0] = Z[J = 0, M = 0; V] = 1. 
Let us define the generating functional W[J, M] — —%hxZ\J,M\, By means of this definition we introduce two 
new functionals: R a (z_; J, M) — SW[J, M]/SJ a (z_) and G(x,y; J, M) = —i5W[J, M]/5M(y,x). When the sources are 
nullified, the new functionals are equal to the average photon field R a (z) and to the average Green's function G{x, y), 
respectively: 

, _ .... SW[J,M]. 1 SZ[J,M}. 5Z[J,M]. 

R a (z; J,M)\j =M =0 — r- t / \ \j=m=q — , \j=M=o r . - \,j=m=o — -i x T , x J=M=0 = -K«U) 

oJ a (z) Z[J,M\ oJ a (z) oJ a (z) 

n, t i\ r\\ 6W[J,M}. 1 5Z[J,M], 5Z[J,M}. 

G{x,yj J,M)\j=m=o - -i sm^U=m=o - z[j^m] ' J=M=0 SM{^xj ' J=M=0 = _ 5M(^) |j = M = = Gfe ^ 

The next step is to derive the SD equations for the corresponding average quantities using the fact that the functional 
measure in (|48|) is invariant under the translations tJj + —> ip + + Sip + , A — > v4 + <5j4. This assumption yields highly 
nontrivial relations among generating functionals and their derivatives which are as follows: 

= J\\(x)(z) + £f|^(z,l')%X)(l'; J, M) - iT ( ^ ±) (y,x | z)G(x,y; J, M), (49) 

G~% x; J, M) = G^~%x) ~ M(y,x) - E(y,x; J, M). (50) 

Here, E(y,x; J, M) is a functional of the sources, but when the sources arc nullified we obtain the mass operator 
T,(y,x) for the average single-particle Green's function: 



E(y,a:) = -r[%,x' | z)i?,|(z) - if ^ | z)R x (z) ~ r\°\y,x' \ z) (~* U=m=o) G-\y^,x) 

- rf\y_,x' | z) (-i dG{ f^ M) h=M=o) G-\y_;,x) - tAK ( ^ .„ ^ ) G ^,(r",fV,*' ). 



(51) 



The term proportional to Ru (z) does not need to be taken into account because of the charge neutrality of the system. 
The next step is to find the relationship between the mass operator and the average two-particle Green's function. 
By solving the SD equations (|49| and (|50|) . one can obtain the sources J a and M as functionals of R a and G. By 
means of the identity: 

_ SJ a (z) _ SJ a (z) 5G(x',j/) t SJ a (z) 8R p (z) 



5M{y,x) 8G[x!,yf) 6M(y,x) 5R p {z) 5M(y,x)' 
we calculate 

6G(x',y;;J,M) i 6G{x> ,y>;J,M) ^ (0) (0) , 



(52) 



In the absence of a symmetry-breaking disorder the —i6G/5M is noting but the two-particle Green's function K. 
However, in the presence of a symmetry-breaking disorder the average generating functional Z[J, M] is not enough to 
obtain the average photon Green's function and the average two-particle Green's function, because they both can be 
written as functional derivatives of Z[J, M] plus terms not directly related to the functional Z[J, M] or its derivatives: 



K 



(53) 



6G(x,y;J,M) 
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D ±iln (z,z') = i 5 . ,J j=m=o + iR±(\\)U; V, J, M)R ±m (z>; V, J, M)\ J=M =o 

= - 5R ^-'f\ M) \j=M=o - iRm\\)(z)Rx(\\){z') + iRxwte V, J, M)R ±m (z'; V, J, M)\ J=M=0 . 
oJ±(\\){z) 

Here, we have introduced the functionals: 

_ _ _ n f\ 6W\J,M;V] „ , _ _ _ ~\ f\ SW[J,M;V] 
G(x,y;V,J,M) = -i 1 ' J , % x) z; V, J, M = 1 ' J 55 
6M{y,x) 0J\\(±){z) 



which depend on the random potential V. Since the average of the product G.G (or R.R) is not the product of the 
averages G.G (or R.R), the term G.G in ([53]) does not cancel G.G. We have already mentioned that by performing 
the disorder averaging of both sides of Eq. we can obtain a relationship between the average photon and the 

average two-particle Green's functions. In other words, the average Green's functions ([55)) and (|5^|) must satisfy the 
equation D ± = L> ( L 0) + D^vf 1 KT^ . Obviously, from the SD equations (gU) and (50]) follows that the terms 
G.G, R.R, G.G and R.R in both sides of the last equation cancel each others. Thus, the cavity polaritons in the 
presence of a symmetry-breaking disorder manifest themselves as common poles of the term — i~jffQr^^-\ J=m=o m 

§M and - ^jfffifl | j=m=o in HMD • 

Next, we rewrite the mass operator £ in the following form: 

E(y,x) = -Tf(y,x' | z)R ± (z; V, J, M)G(x>, y>; V, J, M)\j =M = G- 1 { } /,x) 

- r\°\y,x! | z)K ( J,' I',! ) r(°'(y" j£ " | ^"W^te'.aO 
r?W I *)* f |: ) rf (f,x" I ^W^te'.aO 



The first term in ()56|) is the Hartree term in the presence of a disorder. The second and the third terms link the 
mass operator to the average two-particle Green's function, and therefore, they are the Fock contributions to the mass 
operator. By introducing the vertex function ry(_n : 

K(0) ( V i ) r "W ( ^' 1 ^%X)te',2) = K (yt) r M±)V>^ I ^)D^ x) {z',z), (57) 
where the free two-particle propagator is given by: 

^ (0) (f I') = ^ G( -'- )G( -'-' (58) 
one can rewrite the mass operator in the following form: 

X(y,ZL) = -ri 0) (y,x' | z)R ± (z; V, J, M)G(x>, y'; V, J, M)| J=M=0 G- 1 ( I ', a;) 

+ *r[ | 0) (y,a;' | i)G^)T^,x\ z')D\\(z, z f ) + iT™ (y, x' \ z')G(x', y^)T ± (^,x \ z!)D x (z,z') (5Q) 

4-1 A I 



By comparing the expression for the mass operator (|32[) 
in the absence of a symmetry-breaking disorder with Eq. 
(f56|) . we find two differences. The first one is the pres- 
ence of an additional term, Ea = —ikKG^ 1 . The sec- 
ond difference is related to the corresponding Hartree 



terms. In the presence of a symmetry-breaking disorder 
the Hartree term depends on th e average of the product 
R±(z; V, J, M)G(x', y'; V, J, M) of two random function- 
als, defined by Eq. (|5B"j). 

Let us for a moment replace the average R±G with the 
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product of the averages R±.G: 

R ± (z; V, J, M)G{x, m V, J, M) | j =M =o -> i?±(l)G(x,y). 

(60) 

Although the average of the product of two random quan- 
tities is not equal to the product of their averages, the re- 
placement (|60l) greatly simplifies the equations and allows 
us to map them to the Zittartz's equations. Note, that in 
his work Zittartz took into account only the lowest-order 
contribution from the disorder to the mass operator £a, 
which corresponds to the replacement of K by the free 
two-particle propagator K^°\ In this approximation we 
calculate for the mass operator: 

x) = -rf (y, x\z)R ± (z) + S(t - r')AG fj (r, t', r, t) 



+ *r ( °\y,x> | z')G(£',y')r||(y',£ I z')D l{ (z,z'). 



(61) 



The first and the third terms in (|6"Tj) represent the pho- 
tonic and excitonic contributions to the order parameter. 
By nullifying the sources in Eq. (f^|) we obtain a rela- 
tionship between the average photonic field R± and the 
average single-particle Green's function: 



R ± (z) = iD^(z,z')T^>(y,x \ d)G{x,y) 



(o), 



(62) 



The last equation leads to an equation, similar to Eq. 
(HZ}: 

i? ± (k) = 5>(k-q) J ^G CT (q;o;)/(0 2 (k)-/i 2 ). 

The sum of photonic and excitonic contributions to the 
order parameter is: 



A ( k ) = E 



V(k - q) 



(64) 

Since our random potential has a variance (|33[) , the equa- 
tion (21) of Zittartz assumes the following form: 



(65) 



The Fourier transform of the single-particle Green's func- 
tion has the form (Eq. (24) of Zittartz): 



V G cv [k,uj) G vv (k,uj) 



1 

= -D X 



lo - G vv {lj) +ji(k) + t0+ A(k) + G cv (oj) 

A{k) + G cv {lu) to - G cc {uj) - 77(k) - z0 H 



(66) 



where 



,D = [A(k) + G cv {lo)] 2 - (w - G cc (co) - rj(k) + 4 0+)x 
(w-G l ,.„H+7 ? (k)-z0+). 



(67) 



Here, 7j(k) depends on the chemical potential, and is 
defined self-consistently by the solution of the following 
equation: 



r?(k) = 



l^y(k-q) f g[G cc (q,c) + G TO (q,( 



(68) 



Equations (|M)) - ([6^|) form a closed set of equations, 
that can be solved at any density. In the high density 
regime we can neglect the small corrections to rj due 
to G cc and G vv , and map our equations (f64|) - ([68|) to 
the zero-temperature version of Eq. |T]). As a result 
one could come up with the conclusion that the order 
parameter and the energy gap are gradually suppressed 
up to a critical disorder strength. 

Strictly speaking, RG ^ RG, and therefore, all results 
obtained by using approximation (|60p should be consid- 
ered questionable. More importantly, the RG term does 
not allow us to prove the existence of the Goldstonc 
mode below the critical temperature, as we did in Sec II. 
Going beyond the assumption (|60p is a very challenging 
task, which requires to take into account diagrammati- 
cally irreducible vertex parts and an infinite number of 
diagrams neglected by the assumption RG — RG. 



IV. SUMMARY 



We have applied the CPT Green's function formalism 
to the problem of cavity polaritons in the presence of a 
symmetry-breaking disorder. In contrast with the non- 
linear sigma-model and the saddle-point equations, the 
Hartree term in the mass operator does not allow us to 
map the SD equations to the corresponding equations in 
the work by Zittartz. 

The saddle-point approximation and the replica trick 
could be responsible for the different expressions for the 
mass operator in the presence of a disorder. In the ab- 
sence of a symmetry-breaking disorder, the saddle-point 
approximation leads not only to the correct gap equa- 
tion in the high density regime, but by investigating the 
Gaussian fluctuations about the saddle point one can ob- 
tain the collective mode spectru m 10 ' 19 as well. Thus, we 
might suggest that the replica trick is not the right tool to 
perform the averaging over the random potential in the 
case of cavity polaritons in the presence of a symmetry- 
breaking disorder. 
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